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In this paper, we study some quantum properties of a superposition of displaced squeezed 
two-mode vacuum and single-photon states, such as the second-order correlation function, the 
Cauchy-Schwartz inequality, quadrature squeezing, quasiprobability distribution functions 
and purity. This type of states includes two mechanisms, namely, interference in phase 
space and entanglement. We show that these states can exhibit sub-Piossonian statistics, 
squeezing and deviate from the classical Cauchy-Schwartz inequality. Moreover, the amount 
of entanglement in the system can be increased by increasing the squeezing mechanism. In 
the framework of the quasiprobability distribution functions we show that the single-mode 
state can tend to thermal state based on the correlation mechanism. Generation scheme for 
such states is given. 

PACS numbers: 42.50.Dv,42.50.-p 



I. INTRODUCTION 

Developing new states beside the traditional ones is an important topic in quantum optics 
and quantum information theories. Fock state |n) and coherent state \a) are the most commonly 
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used states in these theories. The single-mode squeezed states of electromagnetic field are purely 
quantum states since they have less uncertainty in one quadrature than the vacuum noise level 



HI 



Additionally, these states exhibit a variety of nonclassical effects, e.g. sub-Poissonian statistics 

n 

and oscillatory behavior in the photon- number distribution [2] . These states can be generated via a 
degenerate parametric amplifier . The third type of states given in the literature is the two- mode 



squeezed states 



4[, which contain quantum correlations between two different modes of the field. 



The importance of these states comes from their connection to the two-photon nonlinear processes, 
e.g. non-degenerate parametric amplifier 5j. These states have been used in the continuous- 
variable teleportation [6], quantum key distribution [jj], verification of EPR correlations {sj], etc. 
The single-mode state-obtained from the two-mode squeezed state by tracing out the other mode- 
cannot exhibit squeezing Q, • Precisely, because of the correlation between modes in the two- 
mode squeezed operator, the squeezing of the quantum fluctuations does not occur in the individual 
modes but it occurs in the superposition of the two modes. 

A great attention has been devoted to produce mesoscopic superposition states. These states 
have interesting distinct characteristics than the classical ones such as interference in phase space, 
squeezing and quantum entanglement [ll]. These remarkable properties present the mesoscopic 
superposition states as powerful tools in quantum information processing, metrology [12I and 
experimental studies of decoherence [jjj]. The most famous superimposed state in the literature is 
the Schrodinger-cat state [la]. There are several proposals for generating superposition of optical 



131 ] and the references 



coherent states in the literature. For recent review the reader can consult 
cited therein. Besides the Schrodinger-cat states various types of superposition have been de- 
veloped, e.g. the superposition of squeezed and displaced number states without [lfj and with 



thermal noise 



171 ]. Moreover, the superposition of multiple mesoscopic states is given in 



and has been generated using resonant interaction between atoms and the field in a high quality 
cavity. The superposition of the two-mode states is discussed in [lj]]. These states-under certain 



1-states [11[ and they can be generated by a 



201 ] . The entanglement of a superposition of 



conditions become very close to the well-known Be 
resonant bichromatic excitation of N trapped ions 
two bipartite states in terms of the correlation of the two states constituting the superposition has 
been discussed in [211 ] . 

Developing new states is an important topic for understanding the boundary between the classi- 
cal and quantum mechanics as well as to cover the needs of the progress in the quantum information 
theory. Moreover, the investigation of nonclassical effects of the quantum states is of considerable 
and continuing interest, since it plays an important role both fundamentally and practically in 
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the quantum information theory. Throughout this paper we study the quantum properties of the 
superposition of squeezed displaced two-mode number states (STDSN), in particular, the vacuum 
and single-photon states. In these states the squeezing mechanism is involved via non-degenerate 
squeezed operator. These states are different from the superposition of the single-mode states 
[la. \v\ in the following sense: They include two mechanisms: (i) entanglement and/or correlation 
between the two-modes, (ii) Two-mode interferences in phase space. These states can be generated 



via two-mode trapped ions 



221 ]. as we will show in section VI. For STDSN we study the single-mode 



second-order correlation, Cauchy- Schwartz inequality, quadrature squeezing, quasiprobability func- 
tions and purity. We show that the nonclassical effects are remarkable in the different quantities. 
Also the single-mode state tends to the thermal state based on the correlation mechanism and the 
amount of entanglement can be increased by increasing the squeezing mechanism. 

We perform this investigation in the following order. In section [U we introduce the state 
formalism and comment on its photon-number distribution. In section IIIII we discuss the second- 
order correlation function and Cauchy-Schwartz inequality. In section HVl the quadrature squeezing 
in the framework of principal squeezing is investigated. In section |V] quasiprobability distribution 
functions and the purity are investigated. The generation of the STDSN is discussed in section IVT1 
however, the conclusions are summarized in section IVlII 



II. STATE FORMALISM 

The correlated two-mode squeezed states are connected with the two-mode squeeze operator, 
which has the form: 

S(r) = exp[ T -(ab-an% (1) 

where a (a') and b (&') denote the annihilation (creation) operators of the first (signal) and 
second (idler) mode, respectively. By means of this operator and the superposition principle we 
develop a new class of states, namely, superposition of squeezed displaced two-mode number states 
(STDSN) as: 

\ip) = \r,a,P) e = \ t [D(ai,a 2 ) + eD(-a 1 ,-a 2 )]S(r)\n,m), (2) 

where e = |e| exp(i(fi), S(r) is given by ([I]) and D(ai, a 2 ) is the two-mode displaced operator defined 
as: 

D(ai, a 2 ) = D{a\)D{a 2 ) = exp(ai<r — a*a) exp(a 2 tr — a 2 b) (3) 
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and ay is generally a complex parameter (a field amplitude), however, throughout the investigation 
in this paper it will be considered real. Also the prefactor X e is the normalization constant, which 
can be easily evaluated as: 

| A e | " 2 = 1 + | e | 2 + 2 1 e | ^L n (4i 2 ) L m (4*1) cos 0, (4) 

where 

t\ = a\C r + a-zS ri t 2 = oliC t + a\S r , S r = sinhr, C r = coshr, \i = exp[— 2{t\ + (5) 

and L n (.) is the Laguerre polynomial of order n (see (8) below). Throughout the paper we study 
only two choices for the parameter [e|, namely, 1 and 0, however, for the parameter <p we take the 
values 0,7r and n/2. Precisely, when |e| = 1 and <fi = 0, 7r,7r/2 the states fl2J) are called even-type, 
odd-type and Yurke-type states, respectively. 

When \n, m) = |0, 0) the states ([2]) can be expressed in a closed form in terms of the Fock states 



2J as: 



^ C(ni,n 2 )|ni,n 2 ), (6) 

ni,n 2 =0 



where 



C(m,n 2 ) = A £ [l + (-ir^ 2e] _2_ eM _l {ai[li +a2/i2)]v /M! (/ii) „ 1 -A/ (/i2) n 2 -M 

x(tanhr)Xr- M (-S), ( 7 ) 

//i = «i — a 2 tanh r, /i 2 = &"2 — c^i tanh r, M = min(ni, n 2 ), TV = maar(ni, n 2 ) 
and Ljj(.) is the associated Laguerre polynomial having the form: 

W-L^!),^!),!,- ( § ) 
The photon-number distribution of © can be evaluated as: 

P(m 1 ,m 2 ) = \C(m 1 ,m 2 )\ 2 , (9) 

where C(mi,m 2 ) is given by ([7]). It is obvious that P(mi,m 2 ) can exhibit pairwise oscillations 
based on the values of the sum mi + m 2 , even if r = 0. We have to remark that the components 
of the STDSN can exhibit oscillatory behavior in P(mi,m 2 ) 23|, apart from the superposition 
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mechanism, which can make this behavior more or less pronounced. Moreover, the single-mode 
photon-number distribution can be obtained via the relation: 

oo 

P(mi) = \C(m 1 ,m 2 )\ 2 . (10) 

rri2=0 

In P(mi) the occurrence of the oscillatory behavior results from the interference mechanism. We 
can explain this fact for the simplest case r = 0,n = m = 0,e = exp(i(j>) and hence (jlpp reduces to 

2mi 

P(mi) = 2|A e | 2 exp(-a 2 )-^-p[l + (-l) mi exp(-2a^) cos <j>}. (11) 

The oscillatory behavior in P(m\) depends on the values of a 2 and (ft, i.e. for large values of 
P(mi) tends to that of the coherent state. This means that one can use the second mode to control 
the nonclassical effects in the first mode and vice versa. 

In the following sections we investigate the properties of the state ([2]) . For the sake of simplicity 
we treat the second-order correlation function, Cauchy- Schwartz inequality, squeezing and purity 
using the form © (,i.e., |n, m) = 1 0, 0) ) , however, the quasiprobability functions are given for the 
case \n,m) = |0, 1). This is to estimate a global information on the generic form. 



III. SECOND-ORDER CORRELATION FUNCTION AND CAUCHY-SCHWARTZ 

INEQUALITY 

In this section we investigate the behavior of the second-order correlation function and Cauchy- 
Schwartz inequality for the state ([6j). These two quantities can give information on the correlation 
between the modes in the quantum system. The second-order correlation function for the first 
mode, e.g. a, is defined by 



where g^(0) = for Poissonian statistics (standard case), g^ 2 \0) < for sub-Poissonian statistics 
(nonclassical effects) and g^(0) > for super-Poissonian statistics (classical effects). The second- 
order correlation function can be measured by a set of two detectors |24|], e.g. the standard 
Hanbury Brown- Twiss coincidence arrangement. For this quantity we restrict the discussion to the 
first-mode only. For this mode one can easily obtain: 



(a f a) = |A £ | 2 {(S 2 + a?)(l + |e| 2 ) + 2|e| / ucos <p [S 2 - 4txt 2 S r C r + 2 ai t 2 S r - 2a x txC r + a\\ }, (13) 
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FIG. 1: The second-order correlation function of the first mode for (|e|,</>) = (l,7r) against {011,0.2) ( a )~(b) 
and r (c) for different values of the parameters as indicated. 



{a) 2 a 2 ) = \K\ 2 {(a\ + 2S? + 4a 2 S 2 )(l + |e| 2 ) + 2|e|Atcos <f> [\a\C 2 T t\ + lQS 2 C 2 t\t 2 2 
+\Qa 1 S r C 2 t 2 l t 2 + a\ + ka\S 2 t\ + 2Sf + 8c*i S?t 2 + 4a? S r t 2 - 8ajS 2 t 2 - AafCrh (14) 
-W ai S 2 C r ht 2 - 165 r 3 C r tit 2 - 16of5 P C r tit2 - 8aiS 2 C r ti]}. 

t is worth mentioning that the most general cases for equations (I13p and (j!4H have been given in 
25 ] for the multimode squeezed cat states but with different parameterizations. 
Substituting ([15]) . ([14"]) into ([T2]) and taking r = we obtain 



y n + 1,12 _ 



exp[— 2(af + al)] cos 4>] 2 
[1 + |e| 2 - 2|e| exp[-2(af + a|)] cos 



(15) 



From ()15p it is obvious that the sub-Poissonian statistics can occur only for <\> = it and 2{a\+a 2 ) 
being small. This means that the odd-type state can exhibit nonclassical effects in the framework 
of#( 2 )(0). In this case, the mode under consideration reduces to the standard odd-coherent state 
with the components | ± \/ a 2 + a 2 ,}. The obvious remark is: when the mode a is prepared in the 
vacuum state |0), its g {2) (0) can exhibit sub-Poissonian statistics based on the values of a 2 of the 
second mode. Similar argument can be given to the second mode. This reflects the role of the 
correlation between the modes in the system, which leads to the possibilities of controlling one mode 
by the other one. Now we draw the attention to the general case when the squeezing mechanism 
is involved. We have noted that the even-type and the Yurke-type states cannot exhibit sub- 
Poissonian statistics. Information about g^(0) of the odd-type states is depicted in Figs. l(a)-(c) 
for given values of the parameters a\ and a 2 . From Fig. 1(a) one can observe the occurrence of the 
sub-Poissonian statistics, in particular, for small values of a's. When the squeezing mechanism is 
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FIG. 2: The factor V against the parameter r for (|e|, 4>, a\) — (1, tt/2, 0.5) (a) and (1, n, 0.5) (b) for «2 = 0.2 
(solid), 0.3 (short-dashed curve), 0.4 (long-dashed curve) and 0.5 (star-centered curve). 

involved, the amounts of the nonclassical effects in g^ (0) decrease and eventually vanish for large 
value of r (see Fig. 1(b)). For ot\ = 02 = and state flSJ) reduces to the two- mode squeezed 

vacuum state. In this case we have g^(0) = 1, which is independent of r. This is remarkable 
in Fig. 1(b), which shows sub-Poissonian statistics only when one or both of otj > 0. Fig. 1(c) 
gives the range of the parameter r (for certain values of ctj ) for which the sub-Piossonian statistics 
occur. It is obvious that the smaller the values of ctj, the greater this range is. 

We conclude this section by investigating the intermodal correlations in terms of the deviation 
:'rom the classical Cauchy-Schwarz inequality. Classically, Cauchy-Schwarz inequality has the form 



{hh)<{h)(h}, 



(16) 



where Ij,j = 1, 2 are classical intensities of light measured by different detectors in a double-beam 
experiment. In quantum theory, the deviation from this classical inequality can be represented as 
V < 0, where the factor V takes the form 271 ] : 



V 



1. 



(17) 



(at&ta&) 

Occurrence of negative values in V means that the intermodal correlation is larger than the cor- 
relation between photons in the same mode and this indicates strong deviation from the 
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29(]. 



classical Cauchy-Schwartz inequality. The origin in this deviation is that in the quantum me- 
chanical treatment we involve pseudodistributions instead of the true ones. This implies that the 
Glauber-Sudarshan P- function possesses strong quantum properties Moreover, the deviation 
from the Cauchy-Schwarz inequality can be observed in a two-photon interference experiment 
For completeness, the expectation value (a^Wab) for the state Q can be easily evaluated as: 



(atfiStfe) = |A e | 2 {[S? + (S r C r - ai a 2 ) 2 + (a? + a 2 2 )S 2 ]{l + |e| 2 ) 
+2|e|/icos (f) [S' 2 cosh(2r) — a\a 2 sinh(2r) + (a 2 . + a 2 )S 2 

-2tit 2 [sinh(2r)(2S' 2 + a? + a 2 . + cosh(2r)) - 2a x a 2 cosh(2r)] + 4sinh 2 (2r)i 2 i 2 . 
+2t2[— S r sinh(2r)a2 + a\S r cosh(2r) — a\a.2(ci\C r — a2S r )] 
— 2ti[{a\C T — a2S r )S 2 + {a\S T — a2C r )(S r C r — 0102)] 



(18) 



- 4sinh(2r)(aiC r - a 2 S r )(t 2 - - 2a x a 2 sinh(2r)(tf + t\) + aja?,} . 

The expectation value {b' 2 b ) can be obtained from (|14|) using the interchange a\ < — > a.2- One 
can easily find V = for r = 0. Generally, we have noted that V < only when ay are small 
(see Figs. 2). Fig. 2(a) is given for the Yurke-type state which is identical to that of the two- 
mode squeezed displaced states (cf. (fT4"|) and (fl~8j) for <f> = vr/2). From these figures the strongest 
deviation from the classical inequality occurs for ay (7^ 0) and r small, i.e. the photons are more 
strongly correlated than it is classically possible, and then the curve monotonically increases as r 
increases. As the values of a 2 increase the negative values in the factor V decrease and eventually 
disappear (compare different curves in these figures). Also a comparison between Fig. 2(a) and 
Fig. 2(b) shows that the nonclassical effects occurred in the factor V for the odd- type state are 
greater than those in the state with Yurke-type state. 

IV. QUADRATURE SQUEEZING 

In this section we discuss the quadrature squeezing for the state under consideration. As it is 
well known, the quadrature squeezing can be measured by a homod yne detector in which the signal 



is superimposed on a strong coherent beam of the local oscillator 



Here we use the notion of 
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FIG. 3: Squeezing factor S for (|e|,</>) = (1,0) with (r, a,\,a2) — (0,a,a) (a) and a\ = a 2 = 0.6 (b). 



the principal squeezing 3l|, which can give one form for the single-mode and two-mode cases. In 
this respect, we define the two quadratures in the following forms: 



X = X\ cos v + X<i sin u, Y = Y\ cos v + Yi sin v, 



(19) 



where the subscripts 1 and 2 stand for the first and second mode, respectively, and v is a rotation 
angle. When v = 0,7r/2,7r/4 the quadratures (Qj)J) yield those of the first mode, second mode and 
compound modes, respectively. For the first mode the quadrature operators can be defined as: 



a + 



Yi = — | a - a 



(20) 



Similar definition can be quoted for the second mode via the interchange a — ► b. The quadratures 
(|19p satisfy the following commutation rule: 



X,Y 



(21) 



Therefore, the squeezing factors associated with the X and Y can be expressed as: 
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(22) 



F = ((A X) 2 ) - 1 

= Fi cos 2 v + F 2 sin 2 z/ + F c sin(2z/), 

5 = ((Ay) 2 )-i 

= S\ cos 2 z/ + S 2 sin 2 v + S c sm(2u). 
where Fj,Sj,j = 1,2, c, take the forms: 

Fi = 2S 2 + 4a 2 |A £ | 2 [l + |e| 2 - |A £ | 2 (1 - |e| 2 ) 2 ], 

S 1 = 2S 2 - 4|A e | 2 |e|/i(af + 4tit 2 S r C r ) cos <j> 

-4|A e | 2 [e[^(2tiC r - ai) 2 (cos 4> + 4|A e | 2 |e| 2 // sin 2 (f>) , 

F c = <±a.\a 2 — 2S r C r — 8aia 2 \e\n\X e \ 2 cos (ft — 4|A e | 2 (l — |e| 2 ) 2 «ia2, 



S c = 2C r S r - 8|e||A e | 2 At(2tiC r - ai)(ii5 P + t 2 C r ) cos^ - 16|e| 2 |A £ | V(«i - 2txC r )(a 2 - 2t 2 C r ) sin 2 

The expressions for F 2 ,S 2 can be obtained from F±,Si via the interchange (ai,t\) < — > (ce 2 ,t 2 ). 
The system is said to be squeezed in x-quadrature or y-quadrature if F < or S < 0, respectively. 
When [e[ = or a\ = a 2 = the expressions (f22j) reduce to 



F = 2S r [S r - C r sin(2z/)] , S = 2S r [S r + C r sin(2i/)] . (24) 

It is evident that for z/ = or 7r, we obtain F = S = 2S 2 , i.e. single-mode squeezing does not 
occur for these cases [3^ ]. On the other hand, for z/ = 7r/4 (, i.e. the compound squeezing factor) 
we have: 

F = exp(-2r) - 1, S = exp(2r) - 1. (25) 



Trivial remark from (|25p . squeezing occurs in the x-component only [32|. This is related to strong 
correlation between the two modes. This behavior is reversed for the superposition state (, i.e., 
|e| 7^ 0), where squeezing can occur in the y-component for certain values of at's and r. For 
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(|e|,r) = (1,0) the expressions ([22]) reduce to: 

F = 8|A e | 2 [ai cos v + «2 sinz/] 2 , 

(26) 

S = — 8\\ e \ 2 /j,(ai cos v + 02 sin i^) 2 (cos + 2|A e | 2 // sin 2 </>). 

It is evident that squeezing can occur only in the y-component for the cases (|e|, (ft) = (1, 0), (1, vr/2). 
Moreover, to obtain squeezing from the single-mode case, the mode under consideration should be 
prepared in a state different from vacuum. When the values of ay increase, i.e. the correlation 
between the two modes starts to play a role, the coefficient /i goes rapidly to zero decreasing 
squeezing inherited in the system. This shows how one can control the behavior of one of the 
modes by the other one. In Figs. 3 we give information on the even-type states for the given 
values of the parameters. In these figures we consider < v < tt/2, which is sufficient to obtain 
full information about the different quadratures. Additionally, the behavior of the quadratures 
in the range ir/2 < v < ir is just a mirror image to that in < v < tt/2. In figure (a) we take 
ol\ = ot2 = a and r = 0. From this figure-regardless of the values of ^-squeezing occurs within the 
range < a < 1.5 otherwise S > 0. Furthermore, the minimum value in S is observable around 
v = 7r/4 and a ~ 0.6. On the other hand, we have noted that squeezing mechanism decreases the 
amount of squeezing involved in the system. This obvious in Fig. 3(b), which shows the range 
of r over which squeezing is available, i.e. < r < 0.35. In other words, r = 0.35 is the critical 
value for a± = c*2 = 0.6. This critical value is a dependent, however, we have found difficultly to 
obtain an analytical form for it. Comparison between Figs. 3(a) and (b) shows that involving the 
two mechanisms (i.e., squeezing and superposition) in the system destroys the nonclassical effects 
contributed by each one independently. 



V. QUASIPROB ABILITY DISTRIBUTION FUNCTION 

Quasiprobability distribution functions, namely, Husimi function (Q), Wigner function (W), 
and Glauber P functions [33:], are important tools in quantum optics. Knowing these functions, all 
nonclassical effects can be predicted and the different moments of the operators can be evaluated. 
Most important, these functions can be measured by various means, e.g. photon counting experi- 



ments 341 ]. using simple experiments similar to that used in the cavity QED and ion traps 35l. l36l]. 
and homodyne tomography j^. In this section, we investigate the single- mode quasiprobability 
distribution functions, in particular, W and Q functions as well as the purity. We start with the 
symmetric characteristic function C w {(3) of the first mode, which is defined as: 
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C W (J3) = Ti\pexp(l3tf -fa)], 



(27) 



where p is the density matrix of the system under consideration. It is mentioning worth that the 
moments of the bosonic operators in symmetric form can be evaluated from C w {(5) by differentia- 
tion. From (J2j) and (|27p one can easily obtain: 



C w ((3) = |A e | 2 |{exp (-^-cosh(2r)) [exp((/3 - /3*)ai) + |e| 2 exp((/3* - $)<*{)] L m (S 2 \(3\ 2 )L n (C, 



(28) 



where 



exp(-i0 - fc+ ^ fc+ )L m (fc / + )L n (fc + ) + exp(i^) - fc ^ )L m (fc'_)L n (fc_) }, 



fc ± = |/JCr ± 2ti| 2 , k' ± = \j3*Sr ± 2t 2 \ 2 (29) 

and L m (.) is the Lagurre polynomial, which can be obtained from (JSj) by simply setting v = 0. The 
W and Q functions can be evaluated, respectively, through the following relations: 

W{z) = ir- 2 fd 2 (3C w (f3)exp{z(3* - /3z*), 

(30) 

Q(z) = it- 2 J d 2 (3C w (P) exp(z(3* - (iz* - ±|/3| 2 ), 
where z = x+iy. Generally, it is difficult to obtain closed forms for these functions for m/0,B^0, 
however, the integration can be numerically treated. Therefore, we restrict ourselves to the case 
n = 0, m ^ 0, which is sufficient to obtain information on the system. On substituting (|28|) into 
(|30p and carrying out the integration we arrive at: 



W(z) 



2|A e | 



7rcosh m+1 (2r) 



452 



2|^+ai| 2 



452 



~ ' _ cosh(2r)J""U- a)s h(2r)l 2: ~ "iTJ + l £ | Z eX PH cosh(2r) ] L ™ [~ cosh(2r) l Z + ai l 



+ 2 H eX P[-^iR20 ( Q 2 + ^ + y 2 )] Re | eX P ( -*> + ^^27)) L mW I } 



4j/A \ 



(31) 



Q(z) = ^^{exp[-^^]L m [-|z - ai \ 2 tanh* r] + \e\ 2 exp[-^^]L m [-|z + atf tanh^ r] 



|2 +„„ u2 , 



k+«i| 2 



|2 i„„i,2 , 



+2|e|exp[-^(a 2 + 2t 2 + x 2 + y 2 )]Re exp (-%<f>- i^j L m (h') }, 



where 



A — t\C r + i 2 »5r> h 



cosh(2r) 



+ % a 2 sinh(2r) 



2ui2l 



*Sy I £ 



(32) 



/»' = ^[(a 2 C r +t 2 ) 2 + 2iyS r (2a 2 C r + aiS r ) - \z\ 2 S 2 ]. 




FIG. 4: The W function of the STDSN for (a u a 2 , \e\,<j>) = (2,0.9,1,0) with (r, m) = (0,0) (a), (0.4,0) (b), 
(0, 1) (c) and (0.4, 1) (d). 



In the derivation of ([31~]) we have used the generating function of Laguerre polynomial [38( namely: 

^^=f>L n (y), (33) 

and the following identity 391 ]: 



/ exp [-B\p\ 2 + (c/2)/3* 2 + ( Cl /2)(3 2 + 7l /3 + 7/5*] d 2 /3 

(34) 

= ^exp{^[7 7lJ B + 7 2 ( Cl /2) + 7l 2 (c/2)]} , 

where X = B 2 — cc\ if Re[-B + ^(c + ci)] and ReiT > 0. In the Appendix we show that the 
quasiprobability functions (|3ip are normalized. It is mentioning worth that the explicit analytical 
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expressions for the W and Q functions of the even/odd superpositions of two-mode squeezed 
coherent states, which are special cases of (I3ip by simply setting m = 0, were obtained earlier in 
[40! . We start the investigation with the W function. The W function has taken a considerable 



34, 



35, 36, 371, and 



interest in the literature since it can be implemented by various means, e.g., 
it is sensitive to the interference in phase space, as we shall show below. From (I3ip we can extract 
several analytical facts. For instance, when (r, m) = (0, 0) and the value of «2 is very small, the W 
function of the first mode exhibits the well-known shape of the cat-state function, i.e., two Gaussian 
bell and interference fringes in-between (we have checked this fact). This can be easily understood, 
where-in this case-the second mode is very close to the vacuum state and hence ([2]) reduces to 
\ip) ~ A e [|ai) +e| — ati)) |0), i.e. the first mode evolves in the Shrodinger-cat state. Additionally, 
when «2 increases the negative values in W function gradually decreases and eventually vanishes 
showing the W function of the statistical mixture of coherent states. This is related to the fact 
that the interference term in the W function includes the factor exp(— 2a|), which tends to zero 
for large values of «2- On the other hand, when |e| = 0, (r, m) 7^ (0,0), the W function cannot 
exhibit neither negative values nor stretching contour in phase space since L m (— k) > as k > 0. 
In this case, the behavior of the W function is close to that of the thermal state for which the peak 
occurred in the W function is greater than that of the coherent light. 

Now we draw the attention to the general case (see Figs. 4 for even-type state). For m = n = 
we have numerically noted that the W function exhibits negative values only when ol\ > a.2 ■ This 
condition can be analytically realized as follows. From (|3ip the interference term in the W function 
includes cos(0 — ^^^r) ) ! wn i cn i s responsible for the occurrence of the negative values in the W 
function. Assuming that we choose the values of the interaction parameters to verify cos(.) = — 1 
and take (|ej, rc) = (1,0) to simplify the problem. Thus, the W function reduces to 

W(0,y)= 4|A „f -{exphS^^-exphMlA}. (35 ) 
v ,y ' vrcosh(2r) 1 ^ L cosh(2r) J ^ L cosh(2r) u v ' 

The W function involves negative values when W(0, y) < 0. Solving this inequality yields the above 
mentioned condition. From Fig. 4(a) one can observe that the W function has two-Gaussian bell 
and interference fringes in between but with negative values smaller than those of the standard 
cat states, e.g. [41]. These fringes can be amplified for certain values of the squeezing parameter 
(compare Figs. 4(a) and (b)). It is mentioning worth that the amplification of the cat states in 
the parametric down conversion has been discussed in [42|. Now we draw the attention to the case 
in which the second mode includes Fock state |1) (see Fig. 4(c)). From this figure the W function 
exhibits two-Gaussian bell around (x,y) = (±ai,0) and inverted peak in-between with maximum 
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FIG. 5: The W function (a), (b) and the Q function (c), (d) of the STDSN for (ai, a 2 , \e\, (j>, m) = 
(2, 0.9, 1, 7T, 1) with r = 0.4 (a),(c) and 1.8 (b), (d). 

negative value. Comparison between Figs. 4(a) and (c) shows that the existence of Fock state in 
the second mode increases the amounts of the nonclassical effects in the first one. Precisely, the 
interference in phase space in the first mode can be controlled by the information involved in the 
second mode. In this respect the nonclassical effects can be transferred from one of the modes 
to the other through the entanglement process. Furthermore, involving the squeezing mechanism 
in the system smoothes out the negative values in the W function (compare Figs. 4(c) and (d)), 
which vanish for large values of r, we get back to this point shortly. In this case we found that the 
W behaves quite similar as that of the thermal light. This is related to the correlation mechanism 
in the system. Comparison between Figs. 4(b) and (d) shows that the squeezing mechanism 
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makes the interference fringes more or less pronounced based on the value of m. Now we draw 
the attention to Figs. 5 given for W and Q functions, as indicated, for the odd-type states. From 
Figs. 5(a) and (b) W function exhibits negative values and with more structure compared to those 
of the even states (compare these figures with Figs. 4(a) and (b)). From Fig. 5(b) the negative 
values still exist even for large values of r. We use the expression "the large values of r" when 
r > 1. This is inspired by the fact: one of the two-mode quadrature operators exhibits maximum 
squeezing (, i.e., F ~ —1 (c.f. (|25p )) when r > 1. From Fig. 5(c) one can observe that Q function 
exhibits a symmetric two- pea k structure, which is representative to the cat states as well as the 
statistical mixture states [4l| . Moreover, when the value of r increases, i.e. the entanglement 
between the two modes becomes stronger, the Q function exhibits a quite similar shape to that 
of the thermal light (see Fig. 5(d)), which is a single peak localized in the phase space origin 
with contour greater than that of the vacuum state. In the framework of W (Q) function the first 
mode exhibits (nonclassical) super-classical light (compare Figs. 5(b) and (d)). This confirms the 
fact that: the W function is more informative than the Q function. Similar conclusions have been 
noticed for the case (|e|, (f>) = (1, 7r/2). We conclude this part by investigating the relation between 
the occurrence of negative values in the W function and the value of the squeezing parameter r. 
To do so we plot Fig. 6 for the W function in terms of r for the same values of the parameters as 
in Fig. 4(c) (dashed curve) and Fig. 5(b) (solid curve). The values of x and y have been chosen as 
they give maximum negative values in Fig. 4(c) and Fig. 5(b). From Fig. 6 we can obtain a rough 
information about the minimal value of r for which the negative values in W function vanish, 
e.g. for even-type and odd-type states it is r = 1 and r = 2.3, respectively. We have checked 
the behavior of the W function for these values and found that the negative values are negligible. 
Furthermore, after plotting the Q function for various values of r (not detailed here) we observed 
that the exact minimal value for even-type state is r = 1.5. Nevertheless, for the odd-type state 
we found that the negative values-even they are very small-are still observed for all values of r. 

Entanglement is a global property of a system. For a bipartite pure state it has been proved 
that there is a unique measure of the entanglement, which is the von Neumann entropy of the 
reduced state of either of the parties [43]. On the other hand, the purity, which gives information 
on the mixedness in the system, can be used to estimate some information on the entanglement 
in the system. In this respect, we can mention that the purity and the von Neumann entropy can 
give quite similar behavior for the quantum system 44|. Also for the Jaynes-Cummings model it 



has been shown that the von Neumann entropy and purity are equivalent 



45J. As the purity is 



easy to be calculated and can provide some exact information about the system, we use it here to 
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FIG. 6: The W function against r for [a.\, |e|, m) = (2, 0.9, 1, 1) with (x, y, <j>) = (0, 0, 0) (dashed curve) 
and (3, 0, tt) (solid curve). The grid is given to show the bounds of r for which W = 0. 



study the mixedness and/or the entanglement in the state under consideration. The single-mode 
purity can be evaluated via the characteristic function through the relation: 



7T 



\C w (P)\ 2 d 2 f3, 



(36) 



where pj is the density matrix for the mode under consideration. For pure (mixed) state we have 
Trp 2 = l(< 1). From ([28]), J5BJ) and setting (|e|, m,n) = (1,0,0) we obtain: 



Tr P? = rJh(£) i 1 + cos ( 2 ^) + ex P (" 



lof 



cosh(2r) 



(37) 



+Ap cos (j) exp 



(01—02X01— A) 

cosh(2r) 



-^ 2 ex P (^y)}. 



In Figs. 7 we have plotted Trpf for the even-type states. From Fig. 7(a) it is obvious that 
for at\ = or a<i = the two modes are disentangled, where Trp| = 1. When ay increases the 
first mode abruptly tends to the partial mixed state (i.e., Trp^ = 0.5), which indicates strong 
entanglement between the two modes. In this case the behavior is quite similar to that of the 
thermal light with mean-photon number n = 1, which satisfies the inequality < Trp 2 < 1. 
This behavior can be analytically realized by evaluating the limiting case (01,02) = (00,00) for 
the purity ([37]) . which gives: 

1 



Trpj 



2 cosh(2r) 



(38) 
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FIG. 7: The purity of the first mode for (|e|, 0) = (1,0) with r = (a) and 0.5 (b). 



It is evident that Tvp\ = 0.5 for r = 0. When the squeezing mechanism is involved in the system 
the degree of mixedness and/or the amount of entanglement is increased and the purity becomes 
more structured (compare Fig. 7(a) and (b)). The purity tends to steady state for large values of 
az's. The value of the purity at (01,02) = (0,0) can be easily obtained from ([3?]) as: 

Trp? = 1 ■ (39) 
cosh(2r) 

It is evident when the value of r increases the amount of mixedness increases, too. From the 
expressions (|38p and (|39p one can realize that the minimal value of the purity can be achieved for 
large r; e.g. for (a,\, a.2,r) = (2, 2, 5) we have Trpf ~ 0.006. 



VI. STATES GENERATION 



In this section we give a generation scheme for the states (j2J) in the frame work of the trapped 
ions. To do so we consider a two-level ion of mass M moving in a 2 — D harmonic potential of 
frequency u> x in the x-direction and uj y in the y-direction. Also a (a') and b (b') represent the 
annihilation (creation) operators for the vibronic quanta in the x- and y-directions, respectively. 
Then the position operators are given by x = Axq(cl + a*), y = Ayo(& + b ), where Axo = 
(2uj x M)~ 2 , Ayo = {2oj y M)~2 the width of the harmonic ground state. Six beams are used to 
drive the interaction with the ion in the cavity; two are propagating in the x-direction detuned 
by ±u x from the transition frequency of the ion. Two are propagating in the y-direction detuned 
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by ±uj y and two are propagating in the x — y plane detuned by ±(uj x + uj y ). Thus the interaction 
Hamiltonian can be written in the form: 



H int = -(jl.E-a^ + h.c), (40) 

where 

E~ = Ei exp[i((w - oj x )t - k\x + + E 2 exp[i((oj + uj x )t - k 2 x + # 2 )] 

+£3 exp[i((uj - Lo y )t - k 3 y + #3)] + £4 exp[i((w + u y )t - k^y + i? 4 )] (41) 

+£5 exp[i((a;o - Wj. - u^f - fc 5 x - k' 5 y + i? 5 )] + £ 6 exp[z((u;o +u x + v y )t - k 6 x - k' 6 y + i? 6 )], 

with a± are the Puali spin operators, Ej^kj^j are the amplitudes, wave vectors, phases of the 
driving modes, and loq is the ionic transition frequency. Using the operator forms for x, y and 
providing that, the field is resonant with one of the vibronic side - bands, then the ion-field 
interaction can be described by a nonlinear Jaynes-Cummings model [46]. In the interaction 
picture and in the Lamb-Dicke limit, it is sufficient to keep the first few terms. Thus we have the 
following effective Hamiltonian: 

Hint = Hi + H 2 , 

Hi = -{gxtf + g 2 a + g$ + g 4 b)&- + h.c, (42) 
H 2 = -{g<r>ab + g§a)w)<7 + + h.c 

where 



9j = ifyrjj exp(itfj - ^rjj), j = 1, 2, 3, 4 



(43) 



9j = -^jVjVj exp[«?j - + 1]'/)}, j = 5, 6, 

with Qj = fj,Ej and rjj = /c 2 (Ax ) 2 ,j = 1,2,5,6, rfi = fc|(Ay ) 2 ,j = 3,4,5,6 and stand for 
the Rabi frequency and the Lamb-Dicke parameter. The motional and internal dynamics can be 
described in the last Hamiltonian by adding other interactions as discussed in to end up with 
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Hi = -\{9i + 9l)a ] + (g*i + 92)a + (g 3 + gift + (g* 3 + 94)b](&- + &+), 

(44) 

H 2 = -[(55 + 9e)ab + (g* 5 + g e )tftf}(&- + &+). 

Under this Hamiltonian any particle prepared in the state ^(| e ) + \g)) w hl stay in this state and 
the dynamics is reduced to that of the motional degrees of freedom only. Now assuming that the 
system is initially prepared in the following state: 

|*(0)} = (|e} + | 5 })|n,m), (45) 

where |e), \g) denote the excited and the ground state of the ion. Also we have dropped the 
normalization constant in (|45p since it has no effect in the following calculations. It is worth 
mentioning that the Fock state |n) can be prepared with very high efficiency according to the 



recent experiments 



481 ] . We proceed, by applying the Hamiltonian H 2 on the state (|45p for a 



duration time t\ we get 

= exp(-iff 2 Ti)|¥(0)) = S(r)\n,m)(\e) + \g)). (46) 

Then we apply Hi for a duration t 2 to get: 

|tf 2 > = expH^ir 2 )|*i) = D( ai ,a 2 )S{r)\n,m)(\e) + \g)), (47) 

where r = i{g^ + g%)T\, ot\ = —i(gx + <7 2 ) T 2 and a 2 = —i(g3 +_Ra)t 2 . We choose the polarization 
in the quantized field so that it affects the excited state only [48| and apply the Hamiltonian H\ 
for a duration T3 we arrive at 

|*3>=exp(-iffi73)|* 2 ) = [(D(p h p 2 )\e) + D( ai ,a 2 )\g)]S(r)\n,m), (48) 

where (3\ = ot\ — i(g% + <? 2 ) r 3 and (3 2 = a 2 — i(g 3 + gX)T%. After that we apply a carrier pulse of 
Rabi frequency Q.q, whose evolution operator is 

U{t) = cos(O i)(|e)(e| + \g)(g\) - ism(n t)(exp(ie)\e)(g\ + exp(-i9)\g)(e\), (49) 

to the state ^3) to get 

1*4} = t>(T 4 )|^3} = [£>(/3i,^ 2 )cos(fi r4) -iexp(i9)sm(n T4)D(ai,a 2 )]S(r)\n,m)\e) 

(50) 

+ p(ai,Q! 2 )cos(fioT4) — iexp(—i0) sia(Q T^)D(j3i, f3 2 )] S(r)\n,m)\g). 
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Then detecting the particle in either of its states gives the states ([2]), where we can choose f3\ = 
—Oil, = —OL2- Throughout the investigation of this paper we have considered that the parameters 
ay and r are real. This can be achieved in the above equations by simply setting, e.g., ~&j = or 
7T, j = 1, ..,4 and #5 = #6 = tt/2. 

VII. CONCLUSION 

The superposition principle is in the heart of quantum mechanics, which can produce new states 
having nonclassical effects greater than those attributed to the components. In this article we have 
studied the quantum properties for a new class of states, namely, superposition of squeezed dis- 
placed two-mode number states. Particular attention has been given to the two-mode vacuum and 
single-photon states of this class. These states include two mechanisms: interference in phase space 
and entanglement between the two modes of the system. We have studied the second-order cor- 
relation function, the Cauchy-Schwartz inequality, the quadrature squeezing, the quasiprobability 
distribution functions and the purity. We have shown that the system can exhibit sub-Piossonian 
statistics even if the mode under consideration is in the vacuum state. This reflects the role of 
entanglement in the system. The deviation from the classical Cauchy-Schwartz inequality has been 
investigated showing that the photons are more strongly correlated than it is allowed classically. 
For certain values of e the system can exhibit squeezing provided that the values of a\ and a<i 
are small. From the W function it has been shown that the single-mode state resulting from this 
class can behave as a thermal state as a result of the correlation process. Also for m = n = 
the W function of the first mode provides negative values only when a\ > «2- The interference 
in phase space of one of the subsystem can be controlled by the information involved in the other 
subsystem. Additionally, the squeezing mechanism can make the interference fringes more or less 
pronounced. The W function is more informative than the Q function in the description of the 
quantum systems. For the purity it has been shown when the values of a's increase the mode 
under consideration abruptly tends to the partially mixed state. The amount of entanglement in 
the system is increased when the value of r is increased, too. Also we have discussed how this class 
of states can be generated by means of trapped ions and pulses for appropriate durations. 
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Appendix 

In this Appendix we prove that the W and Q functions (]3ip are normalized. Precisely, we would 
like to prove the followings: 



W(z)d 2 z = 1, / Q(z)d 2 z = 1. 



(51) 



To do so we use the generating function technique. We focus the attention on the Q function only, 
where the W function can be similarly treated. Moreover, we evaluate the integration for one of 
the interference terms in the Q function, which we denote I m and has the form: 



J ™ = J eX P[-^2 («1 + 2t 2 + N )] ex P 



Multiply both sides of (|52h by t m , hence sum over index m and use the identity (|33h we obtain: 



m=0 



exp(— i(f>) 
tt(1 - ti)C? 



exp 



(gf + gg) 



cxp 



|z| 2 A(z - z*) + tit' 
C3 tf 1 



d 2 z, 



(53) 



where t' = t/C 2 . Invoking the value of h' from (f3"2"j) into (f53"j) and apply the identity we arrive 
at: 



^ Im 

m=0 



cxp(-^) 

1-t eX P 



(Ql+ 2t 2) , f , (a2Cr+t 2 ) 2 



73" 



+ 



(t'-i)c? 



r [S r (2 a2 C r +q 1 S,.)t'-A(t'-l)l 2 -i 

XKX P\ (t'-l)(t-l)C2 J' 

The exponent in the above equation can be rewritten in terms of the parameter t as: 



oo 

t I-m 

m=0 



cxp(-i^) 

1-t eX P 



(q^+2t|+A 2 ) _ ^Cr+ai-gr+SrA] 2 * 
C*(t-1) 



exp(— i0) exp 



(a 2 +2t 2 +A 2 ) 
E3 



oo 
m=0 



(2o 2 C r +aiS r +5 r A) 2 



(54) 



(55) 



The transition from the first line to the second one has been done by means of the identity (|33p . 
Now the value of the required integral is 



I m = exp(-i(/>)exp 



(oj + 2t| + A 2 



L, 



(2q 2 C r + ai 5 r + g r A) ; 



(56) 
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Through minor treatments one can easily prove: 

(of + 2t\ + A 2 ) _ 2 2 {2a 2 C r + ai S r + S r A) 2 _ 2 

—2 ~ z \Pl~i~ l 2)i — 4r 2- l°'J 

Therefore, the quantity I m takes the form: 

I m = exp(-*0) exp[-2(t 2 + ^)]L m (4^). (58) 

The value of the integration of the second interference term is just the complex conjugate of (I58p . 
Similar procedures lead to the followings: 



7T f \z — ai| 2 nT r i2 i 2 i t2 ^ /" r k + ^lPix r i i2 ,2 

c 2 m +2 J ex Pi ^2 \Lm[-\z-ai\ tanh rjd z = 2m+2 / exp[ ^ JL m l-p+«i| tanh ? 

(59) 

From these results we conclude: 

J Q(z)d 2 z = |A e r 2 |A e | 2 = 1. (60) 

Using procedures similar to those given above one can prove that the W function is normalized 
(we have checked it). 
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